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ABSTRACT 

The most irradiated transiting hot Jupiters are characterized by anomalously inflated radii, some- 
times exceeding Jupiter's size by more than 60%. While different theoretical explanations have been 
applied, none of them provide a universal resolution to this observation, despite significant progress 
in the past years. We refine the photometric transit light curve analysis of 115 Kepler giant planet 
candidates based on public Q0-Q2 photometry. We find that 14% of them are likely false positives, 
based on their secondary eclipse depth. We report on planet radii vs. stellar flux. We find an in- 
crease in planet radii with increased stellar irradiation for the Kepler giant planet candidates, in good 
agreement with existing hot Jupiter systems. We flnd that in the case of modest irradiation received 
from the stellar host, giant planets do not have inflated radii, and appear to have radii independent of 
the host star incident flux. This flnding suggests that the physical mechanisms inflating hot Jupiters 
become ineffective below a given orbit-averaged stellar irradiation level of ~2xl0* erg s~^ cm^'^. 

Subject headings: planetary systems - techniques: photometric 



1. INTRODUCTION 

The onset of giant exoplanet tran sit science, initiated 
by t h e discovery of HP 209458 b (jCharbonneau et all 
120001: iHenrv et al.l I2OOO0 . immediately revealed unex- 
pected anomalous radii for several hot Jupiters . The 
so-called inflated radii became common as new transit- 
ing hot Jupiters were discovered. The anomalous giant 
planet radii were unexpected because in the mass regime 
of giant planets, the mass-radius relationship for giant 
plane ts was expected to be unique (jZapolskv fc Salpeteii 
Il969f ). assuming a given composition. Moreover, the 
compensating effects of electron degeneracy and elec- 
trostatic contribution from the classical ions yield a 
quasi-constant radius around 1 and 7 Mj^p (see, e.g., 
iChabrier et al.l [2009D. Soon a ft er the discovery o f 
51 Peg b (IMavor fc Quelol HOOl . IGuillot et al.1 (119961 ) 
correctly pointed out that strongly irradiated giant plan- 
ets do not follow the same mass-radius relationship than 
the one of isolated objects. 

More than 100 known transiting hot Jupiters later, 
planet radius discrepancies are stil l common (e.g., 
iBaraffe et alllMoHf^tnev et al.|l2010D. Despite numer- 



ous t heoretical studies fsee.e.g.. iFortnev fc NettelmannI 
120101 for a review), no universal mechanism seems to 
fully account for the observed radii anomalies. 

The motivation of the present work is that most of 
the proposed mechanisms to explain inflated radii are 
expected to become less effective as the stellar incident 
flux decreases (e.g.. IFortnev et all 120071 ). We therefore 
study a Kepler subsample of 138 giant planets candi- 
dates to better understand the effect of irradiation across 
a wide range of orbital separation on giant extrasolar 
planets radii. We further note that transiting giant 
planets receiving modest stellar irradiation are particu- 
larly important for the derivation of their internal struc- 
ture and composition, as the interior energy source is 
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expe cted to affect the planetary radius only nominally 
(Mill er fc FortnevI |20T1 . 

The paper is organized as follows. We flrst describe in 
Sect. 2 how the sample of giant planet candidates was 
chosen. Then we present the data analysis that reflned 
the system parameters and provided constraints on plan- 
ethood of the candidates. The results of this analysis are 
then shown in Sect. 3. We flnally discuss the behavior of 
the Kepler giant planet candidates in the radii vs. stellar 
incident flux plane and estimate the Kepler false positive 
rate for this class of objects in Sect. 4. 

2. DATA ANALYSIS 
2.1. Selection oj giant planet candidates 

This study is based on quarters QO, Ql and Q2 Kepler 
data that were pub licly released on 2011 February 1st 
(jBorucki et al.ll20Tltl . In total, the datasets encompass 
136 days of photometric monitoring between May and 
September 2009. 

A list of 1235 Kepler Objects of Interest (KOI) was re- 
leased, unveiling a huge diversity of exoplanet candidates 
in terms of planetary radii and orbital periods. One of 
the key elements of this release is the relati vely low occur- 
rence of Jupiter-sized candidates (see, e.g.. lHoward et al] 
1201 If ). Out of the 1235 planet candidates, only about 
15% have radii above Q.QRjup, a result that supports the 
low frequency of gia nt planets found in RV surveys (e.g., 
Howard et al.ll2010t ). albeit for a different stellar popula- 
tion. 

Our primary select ion criterion is the planetary radius. 
iBorucki et al.l ()201lD announced 165 giant planet candi- 
dates with 6i?0 < Rp < 15i?0 and 19 candidates with 
15i?® < Rp < 22i?0 for a total of 184 objects. We fur- 
ther restricted this sample to keep only "giant" planet 
candidates, deflned here by 8i?© < Rp < 22i?0. This 
step yielded 138 candidates. We removed the 14 KOI 
that have only one transit in Q0-Q2 and 9 other classi- 
fled as "moderate" candidates, likely false positives, that 
exhibit centroid motion or difference of depth between 
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odd and even transits (jBatalha et alJlMol) . This selec- 
tion left us with a sample of 115 giant planet candidates. 

2.2. Method 

For each of the 115 KOI, we retrieved the Q 0-Q2 
raw long-cadence photometry (jJenkins et al.l[20To| ) from 
MAST0. These data include all photometry in the form 
of individual lightcurves. We used the raw photome- 
try in stead of the Kepler-corrected (PDC) (jJenkins et al.l 
|2010[) photometry so that we could identify systematics 
on specific timescales as necessary input to our analysis. 
Moreover, by using the raw data we can assess the ampli- 
tude of correlated noise (from instrumental systematics 
and stellar variability combined) to derive uncertainties 
on stellar and planetary parameters. 

To better characterize the 115 planet candidates, 
we performed individual Markov Chain Monte Carlo 
(MCMC) analysis for each KOI. The aim of this analysis 
was two-fold. The first goal was to remove false positives 
by way of detecting a robust secondary eclipse signal in- 
dicative of a stellar companion instead of a planetary 
companion. The second goal was to derive the stellar 
density from the transit lightcurve in order to further 
derive the stellar parameters and planetary radius. This 
step also required use of the stellar r,. f f values draw n 
from the Kepler Input Catalog (KIC: lBrown et al.ll20TTh . 

We used the implementation of the M CMC algorithm 
presented in [Cillon ct al. (2009, 120101) in order to de- 
rive the stellar and planetary parameters. MCMC is 
a Bayesian inference method based on stochastic sim- 
ulations that samples the posterior probability distribu- 
tions of adjusted parameters for a given model. Our 
MCMC implementation uses the Gibbs sampler and the 
Metropolis-Hastings algorithm to estimate the posterior 
distribution function of all jump parameters. Our nom- 
inal model is based on a star and a transiting planet on 
a Keplerian orbit about their center of mass. 

Input data provided to each MCMC consisted of the 
Q0-Q2 if ep/er photometry and the KIC stellar T^ff value 
relevant to each candidate. Two runs were performed, 
each of them made of two Markov chains of 50,000 steps 
each. The purpose of the first run was to estimate the 
level of correlated noise in each lightcurve and to provide 
the second run with updated error bars on the jump pa- 
rameters. In the second MCMC, the good mixing and 
convergence of the Markov chains w ere assessed using 
the G elman-Rubin statistic criterion (jCelman fc Rubini 
[l99l) . 

We divided the total lightcurve in chunks of du- 
ration of ^24 to 48 hr and fitted for each of them 
the smooth photometric variations due to stellar vari- 
ability or instrumental systematic effects with a time- 
dependent quadratic polynomial. Baseline model co- 
efficients were determined at each step of the MCMC 
for each light curve with the sin gular value decomposi- 
tion method (jPress et al.l 119921 ). The resulting coeffi- 
cients were then used to correct the raw photometric 
lightcurves. 

For each c hunk of data, correl a ted noise was accoun ted 
for following IWinn et all (l2Q08l ) : IGillon et"all (l2010l) . to 
ensure reliable error bars on the fitted parameters. For 
this purpose, we compute a scaling factor based on 

^ http:/ /archive. stsci.edu/kepler/ 



the standard deviation of the binned residuals for each 
lightcurve with different time bins. The error bars are 
then multiplied by this scaling factor. 

The rest of the important inputs for the MCMC is as 
follows. For each quarter, we estimated the degree of 
photometric dilution by using the contami nation factoin 
computed from the KIC crowding matrix (jBrvson et al.l 
|2010[) and was then applied to the transit photometry. 

We assumed a quadratic law for the limb-darkening 
(LD) and used ci = 2ui + U2 and C2 = wi — 2u2 as jump 
parameters, where ui and U2 are the quadratic coeffi- 
ci ents. ui and U2 wer e draw n from the theoretical tables 
of iClaret fc BloemenI ()2011[ ) for the corresponding KIC 
Teff and log g values. 

The MCMC has the following set of jump parameters: 
the planet/star flux ratio, the impact parameter b, the 
transit duration from first to fourth contact, the time 
of minimum light Tq, the orbital period, the occultation 
depth, the two LD combinations ci and C2 and the two 
parameters i/ecosw and y/e sin uj. A uniform prior dis- 
tribution is assumed for all jump parameters but ci and 
C2, for which a normal prior distribution is used, based 
on theoretical tables. 

2.3. False positive assessment via secondary eclipses 

The Kepler giant planet candidates list is not guaran- 
teed against false positives, although a ranking of prelim- 
inary assessment is provided on MAST. Yet it is the false 
positive rate contributed by eclipsing binary stars with 
larger radii than Jupiter's that would contaminate our 
findings. Furthermore, late type M dwarfs could actually 
produce planet-to-star area ratio indistinguishable from 
bona fide giant planets, beca use of their similar ra dii to 
Jupiter-like objects (see, e.g.. iChabrier et al] |2009'). Ob- 
taining radial-velocity measurements at orbital quadra- 
ture for more than 100 giant planet candidate objects 
is unrealistic for the purpose of the present study, given 
the faint V magnitude of the host stars and the num- 
ber of targets of higher priority in the Kepler follow-up 
program. 

We therefore used our MCMC method to search for 
a secondary eclipse whose depth would be indicative 
of a stellar companion instead of a planet. No con- 
straint on the eccentricity was imposed since binaries 
with orbital periods of a few days only and eccentric- 
ity e > 0.2 are not uncomm on (see, e.g., iRucinski et al.l 
[2007t [Maceroni et al.|[2009l ). 

We used the derived occultation depth ^ to com- 
pute the corresponding geometric albedo Ag — 

and the brightness temperature to assess the nature of 
each KOI. We further visually inspected the individual 
folded lightcurves, as discontinuities due to spacecraft 
roll, change of focus, pointing offsets or safe mode events 
could create artifacts in the raw photometry and affect 
the detection of shallow features in the lightcurve. We 
present the results of this analysis, in terms of false pos- 
itive rate, in Sect. 3. 

2.4. Stellar parameters 

The photometric calibration of the Kepler field target 
stars presented in the KIC yields stellar radii uncertain- 

^ Contamination values can be found in the fits files headers. 
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ties of 35% RMS (|Brown et al.ll201lD . Because the stellar 
radius uncertainty translates directly to a planet candi- 
date size, the stellar radius uncertainty is too large for 
any useful constraint on the behavior of planetary radii 
with incident stellar flux and orbital distance. Hence this 
motivated us to derive our own stellar radii by a differ- 
ent method than assuming the KIC stellar radii, which 
yields smaller uncertainties (of ^15% RMS) on the de- 
rived stellar radius (see Sect. 3). 

The method we use employs the empirical calibra- 
tion law presented infTorrcs ct al. (2010). The authors 
show that accurate stellar masses and radii could be de- 
duced from the stellar effective temperature Tg//, surface 
gravity log g and metallicity [Fe/H] derived from spec- 
troscopy. For this purpose they build a calibration law 
based on a large sample of well characterized detached 
binaries. A linear regression algorithm then provides the 
st ellar mass, in function of the spectroscopic parameters. 

lEnoch et al.l ()201C1[) further suggested to use as input 
the stellar density instead of the stellar surface gravity 
log g. The advantage of this approach is that the stellar 
density is well constrained by the transit lightcurve pho- 
tometry (Seager & Mallcn-Ornclas 2003), yielding better 
results than using the surface gravity derived from the 
spectroscopic analysis. 

The empirical calibration implemented in the MCMC 
is therefore a function of T^ff, Pi, and the stellar metallic- 
ity [Fe/H ], which is poorly co nstrained from the KIC pho- 
tometry (! Brown et al.ll20l"TI ). We thus imposed a 0.3 dex 
uncertainty on the stellar metallicity. At each step of the 
MCMC, Pi, (deduced from the jump parameters), Tg// 
and [Fe/H] (drawn from the normal distribution based on 
the KIC value with the error bars quoted above) are used 
as input to the calibration law. The physical parameters 
of the system are then deduced using the resulting stel- 
lar mass. The intrinsic uncertainty of the parameters of 
the calibration relationship is accounted for by randomly 
drawing the para meters values f rom th e normal distribu- 
tion presented in [Torres et al.l (|20100 at each iteration 
of the MCMC. The remainder of the uncertainty on the 
stellar radius is then mostly dominated by the error on 
the KIC T^f f and on the intrisic scatter of the empirical 
relationship. 

This method makes the derivation of the stellar mass 
possible at each step of the MCMC without the need of 
performing a separate analysis based on stellar evolution 
models. 

3. RESULTS 

The main result of this study is that giant planet can- 
didate radii are independent of stellar incident flux be- 
low an incident flux of about 2 x 10* erg s~^cm~^ (Fig- 
ure 1). Although the giant planet radiu s trend was hinted 
with published giant planets alone (jMiller fc FortnevI 
20111 see also Fig. 1) and theoretically expected (e.g., 
Fortnev et al.ll2007l ). inclusion of these new Kepler giant 
planet candidates radii yields a robust trend. 

The objects supporting this result are the Kepler gi- 
ant planet candidates that have no or shallow secondary 
eclipses consistent with their equilibrium temperature at 
2(7 level or less. For comparison we include transiting 
planets not discovered by Kepler which overlap perfectly 
in the Rpianet vs. incident flux plane, but mainly pop- 
ulate the high incident flux regime. The Kepler Q0-Q2 



coverage not only enables almost a doubling of the tran- 
siting giant planet candidates but also extends the cover- 
age out to lower incident fluxes as compared to currently 
known transiting planets. 

Complicating the result is the fact that Kepler planet 
candidates are not vetted as actual planets. We have 
used our MCMC analysis to assess false positives via 
secondary eclipse detection. Indeed, 16 planet candidates 
show strong evidence for deep secondary eclipses suggest- 
ing a 4cr discrepancy or more with their estimated equi- 
librium temperature. Such objects are discarded from 
the study. Finally, 22 planet candidates yield a secondary 
eclipse signature whose origin cannot be secured, the in- 
ferred brightness temperature being consistent with ei- 
ther a planetary or stellar companion. We still choose to 
include those candidates in Figure 1, with distinct sym- 
bols. Additional data will help in tightening the nature 
of those objects. 

We notice that eclipsing binaries with grazing transits 
combined to non-zero orbital eccentricity would not yield 
any secondary eclipse and would therefore be wrongly 
identified as planets in our study. Any such contamina- 
tion should be uniform with the range of incident fluxes 
explored in this study and would not affect the main 
finding of a trend in giant planet radii. 

For most of our planet candidates, there is no infor- 
mation about the orbital eccentricity. Instead of un- 
realistically assuming circular orbits, we assigned each 
candidate an eccentri city value dr a wn fr om the distri- 
bution presented in iWang fc FordI ()2011| ) , as well as a 
random value for the argument of periastron. This ap- 
proach is reasonable since no significant trend seems to 
exist betwe en eccentricity and orb ital period for Kepler 
candidates (jMoorhead et al.ir2011l ). 

To gauge the impact of orbital eccentricity on our re- 
sults, we performed a new MCMC analysis by imposing 
priors on y/esinuj and ^/ecosu! (see Sect. 2.2), based on 
the values drawn for e and uj in the previous step. We 
then computed the orbit-averaged incident fiux for each 
candidate and found an excellent agreement (4% in aver- 
age) with the fiuxes obtained for the circular case. The 
reported trend in the Rp vs. incident fiux plane is there- 
fore robust to the planetary candidates' orbital eccen- 
tricity. We show our results (assuming the eccentricity 
distribution described above) on Fig. 1. 

In summary, out of the 115 KOI, 70 of them exhibit no 
or shallow secondary eclipses consistent with their equi- 
librium temperature. Those objects are therefore con- 
sidered as of planetary origin, whereas 16 are identified 
as stellar companions and 22 are of ambiguous classifica- 
tion. We list the KOIs and their identification in Table 1. 

In the course of the analysis, 7 KOI got their radius 
revised to less than 6 i?® . Those KOI were thus discarded 
from our study, as no longer part of the giant planet 
candidate sample. 

4. DISCUSSION 

4.1. Giant planet radii inflation mechanisms 

KOI classified as planetary candidates in the frame 
of this study yield planetary radii vs. incident flux in 
good agreement with published transiting planet data, 
as shown on Fig. 1. Remarkably, this set of KOI results 
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TABLE 1 



Planetary Planetary (continued) Ambiguous False Positives 



1.01 


398.01 


12.01 


194.01 


2.01 


410.01 


187.01 


197.01 


10.01 


417.01 


189.01 


208.01 


17.01 


418.01 


458.01 


552.01 


18.01 


421.01 


617.01 


609.01 


20.01 


423.01 


728.01 


743.01 


22.01 


425.01 


763.01 


745.01 


94.01 


625.01 


767.01 


779.01 


97.01 


674.01 


772.01 


876.01 


98.01 


686.01 


823.01 


895.01 


100.01 


698.01 


840.01 


1003.01 


127.01 


760.01 


855.01 


1152.01 


128.01 


801.01 


856.01 


1177.01 


135.01 


805.01 


918.01 


1540.01 


138.01 


806.01 


929.01 


1541.01 


183.01 


806.02 


960.01 


1543.01 


186.01 


809.01 


961.02 




188.01 


815.01 


961.03 




190.01 


824.01 


1020.01 




191.01 


846.01 


1285.01 




192.01 


850.01 


1299.01 




193.01 


858.01 


1385.01 




195.01 


871.01 






196.01 


882.01 






199.01 


883.01 






202.01 


889.01 






203.01 


897.01 






205.01 


908.01 






214.01 


913.01 






217.01 


931.01 






254.01 


1089.01 






351.01 


1176.01 






366.01 


1227.01 






368.01 


1391.01 






372.01 


1486.01 






Note. — 


List of KOI used in 


this study. KOI meeting 


our criteria 



Sect. 3) arc sliown in tlie "planetary" column wliile the probable stellar 
companions are shown in the "False positives" column. The intermediate 
class is shown under "ambiguous" . 




10' 10° 10' 

Incident flux [erg s~' cm~^] 

Fig. 1. — Planetary radii as a function of incident flux. Black 
filled circles are KOI ranked as planetary candidates in the frame 
of this work v^hile gray diamonds represent KOI whose origin is 
ambiguous (see Sect. 3). Transiting giant planets previously pub- 
lished, and mostly from ground-based surveys, are shown as red 
triangles. The relevant parameters Rp, Rg, T^ff and a have been 
drawn from http: / / www.inscience.ch/ transits | on August 29, 2011. 



in constant radii of ~ 0.87±0.12i?j|j, similar to Jupiter, 
below an incident flux of ^ 2 x lO^erg s~^cm~^. We re- 
port no inflated giant planet radii below this threshold. 

Several explanations have been invoked to bring or 
maintain heat in the planetary interiors, necessary to 
explain radii anomalies through a larger equilibrium ra- 
dius. An exhaustive d e script i on can be found in, e.g., 
iFortnev fc NettelmannI (pOlol ): iBaraffe et all (|20lol ) but 
none seems to reproduce all planets with inflated radii. 

Tidal energy dissipation in the giant planet in- 
terior is expected to co unteract the cont r action 
(|Bodenheimer et all [200ll ). iBodenheimer et all (|2003l ) 
proposed that an additional companion in the system 
could pump the planet eccentricity that would be dissi- 
pated through tides. Significant f ollow-up on th i s pos- 
sibility emerged in the last years. iLevrard et all ([20091 ) 
have for instance shown that most known transiting plan- 
ets w ere spiraling toward their star due to tidal dissipa- 
tion. I Miller et al.l ()2009() modeled coupled thermal evo- 
lution and tidal effects on giant planets and were able to 
reproduce the radius anomalies for 35 out of 45 planets 
part of their sample, assuming ad hoc initial conditions. 

Electrical current generated through the interaction 
of ionized particles with the planetary magnetic field 
causes a dissipation of ener g y in the planet a ry int erior 
(jBatvgin fc StevensonllMol) . iLaughlin et al.l (|201lD find 
support for this hypothesis from the set of transiting 
planets known in 20f0 but also state that other pro- 
cesses should be contributing to account for the observed 
anomali es, such as the effec ts of heavy elements abun- 
dances (|Batvgin et al.l [201 It) or the internal heating in- 
duced by tidal circularization for eccentric planets. 

Layered convection should occur in atmo- 
spheres ch aracterized by molecular weight gradients 
(iChabrier fc Baraffc 2007). This would deer ease the loss 
of heat and slow down the contraction in the planetary 
interiors. This mechanism is independent of the incident 
stellar flux. 

Our results suggest that the combinations of mecha- 
nisms responsible for the giant planets inflated radii are 
correlated to the strength of the stellar incident flux. 

For the most close-in planetary candidates, tidal ef- 
fects could contribute significantly and add to the sole 
incidence of stellar irradiation. 

As part of this study, we find no support for a process 
that would be solely based on layered convection as no 
cold inflated planets are reported. 

There is observed scatter in the plateau of giant plan- 
ets radii below ~ 2 x 10^ erg s~^cm~^. We speculate 
this scatter could be due to the effects of metallicity. 
Planets enriched with heavy elements yield a more com- 
pact structure, t hus a smaller radius, like HD 149026 b 
(|Sato et al.ll200,^ . 

4.2. False positives rate 

A byproduct of our MCMC analysis is also an impor- 
tant result: an estimated false positive rate for the Kepler 
giant planet candidate listing. According to our analy- 
sis, 14% of the Kepler giant planet candidates studied in 
this sample are eclipsing binaries or background eclips- 

* Jupiter radius itself is 0.977 Rj when using its mean radius of 
69,894 km instead of its equatorial radius of 71,492 km (see Table 1 
of GuiUot 2005) 
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ing binaries. This result is based on identification of 
a large discrepancy between the candidate's equilibrium 
temperature and its measured brightness temperature in 
the Kepler bandpass. We further note that the 14% false 
positive rate number might even be higher in the case of 
grazing events or eccentric orbits, for which the detection 
of a secondary eclipse is not possible. We found no false 
positives among the 10 multi-planetary systems included 
in our sample. 

Additional Kepler data will allow an extension of the 
present study to planets with longer orbits than the can- 
didates presented in this work. The next quarters of 
Kepler photometry will also help in tightening the exact 
fraction of false positives among giant planet candidates 
by improving the characterization of secondary eclipses. 

In summary, this work presents one of the first results 
from the emerging science of exoplanet statistics enabled 
by Kepler's exquisite photometry and large pool of planet 
candidates. With future Kepler data, we expect many 
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